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Abstract 

We investigate the charm quark system using the relativistic heavy quark action on 2+1 flavor 
PACS-CS configurations previously generated on 32 3 x 64 lattice. The dynamical up-down and 
strange quark masses are set to the physical values by using the technique of reweighting to shift 
the quark hopping parameters from the values employed in the configuration generation. At the 
physical point, the lattice spacing equals a -1 = 2.194(10) GeV and the spatial extent L = 2.88(1) 
fm. The charm quark mass is determined by the spin-averaged mass of the IS charmonium state, 
from which we obtain m^ rm (^L = m^—^) = 1.260(1)(6)(35) GeV, where the errors are due to 
our statistics, scale determination and renormalization factor. An additional systematic error from 
the heavy quark is of order a^f(mQa)(aAQCD), which is estimated to be a percent level if the 
factor f(rriQa) analytic in mQa is of order unity. Our results for the charmed and charmed-strange 
meson decay constants are fjj = 226(6)(1)(5) MeV, fn s = 257(2) (1) (5) MeV, again up to the 
heavy quark errors of order ot 2 s f{mQa)(akQCD)- Combined with the CLEO values for the leptonic 
decay widths, these values yield \V cd \ = 0.205(6) (1) (5) (9), (V^l = 1.00(1)(1)(3)(3), where the last 
error is on account of the experimental uncertainty of the decay widths. 
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I. INTRODUCTION 



Precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 
is an indispensable step to establish the validity range of the standard model, and to search 
for new physics at higher energy scales. Lattice QCD has been making steady progress in this 
direction. For the matrix elements such as \ V u d\ and \ V US \ in the first row which involve only 
light quarks, dynamical simulations including up, down and strange quarks have reached 
the point where the relevant pseudoscalar meson decay constants and form factors are being 
determined at subpercent precision. On the other hand, for \V c d\ and \ V CS \ in the second row, 
the precision of lattice QCD calculation for the decay constants and form factors is still at 
5 to 10% level. This is not clearly superior to non-lattice QCD determinations. Indeed, the 
estimate quoted in Particle Data Group (PDG) 2010, \V c d\ = 0.230(11) |l| with an accuracy 
of 5%, is obtained from neutrino and anti-neutrino experiments. 1 Much effort is needed on 
the part of lattice QCD toward a better precision in the charm sector. 

One of the difficulties with the charm quark in lattice QCD simulations at a typical cutoff 
a -1 f=s 2 GeV resides in significant cutoff errors due to the charm quark mass. The heavy 
quark mass correction is rriQa ~ 1, and hence we must control errors to all orders of rriQa 
to achieve a few percent accuracy. The Fermilab action [2| and the relativistic heavy quark 
action 0, 0] have been proposed to meet this goal. In the present work we employ the 
relativistic formalism of Ref. [3[ to explore the charm quark system. 

Another source that prevents precise evaluations in lattice QCD is the error associated 
with chiral extrapolations in the light quark masses. This problem has been increasingly 
alleviated through progress toward simulations with lighter and lighter dynamical quark 
masses and sophisticated application of chiral perturbation theory techniques. The acceler- 
ation of dynamical lattice QCD simulation using multi-time steps for infrared and ultraviolet 
modes j5|, @] has made it possible to run simulations with light up, down and strange quark 
masses around their physical values 0]. In such simulations, uncertainties due to chiral 
extrapolations are drastically reduced. 

In fact we can proceed one more step and reweight [8| dynamical simulations such that 
dynamical quark masses take exactly the physical values. A potential difficulty with dynam- 
ical lattice QCD is a large fluctuation of quark determinant ratios necessary for reweighting. 
We have demonstrated the feasibility of this procedure in Ref. j9j by reweighting a set of 
PACS-CS configurations with m n = 152(6) MeV and m K = 509(2) MeV to those with 
m n = 135(6) MeV and mx = 498(2) MeV. Once the reweighting is successfully made, am- 
biguities associated with chiral extrapolations are completely removed. In the present work 
we employ the reweighting factors and the set of original dynamical configurations employed 
in Ref. [9j. Hence our light quark masses sit at the physical point. 




In this paper we present our work for the charm quark system treated with the relativistic 
heavy quark formalism [3| on the 2+1 dynamical flavor PACS-CS configurations of 32 3 x 64 
lattice generated with the Wilson-clover quark and reweighted to the physical point for up, 
down and strange quark masses. The lattice spacing is estimated as a -1 = 2.194(10) GeV. 
We measure the masses and decay constants of charmonia, charmed mesons and charmed- 
strange mesons. We then calculate the charm quark mass and the CKM matrix elements. 

This paper is organized as follows. Section II explains our method and simulation param- 
eters. Section III describes our results for the charmonium spectrum and the charm quark 



1 jV^sj is hard to be estimated from neutrino and anti-neutrino experiments, |V CS | = 0.94+Q2g ± 0.13 [l| 
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mass. In Sec. IV, we show our charmed meson and charmed-strange meson spectrum. Sec- 
tion V is devoted to present our pseudoscalar decay constants and the CKM matrix elements. 
Our conclusions are given in Sec. VI. 



II. SET UP 

Our calculation is based on a set of Nf = 2 + 1 flavor dynamical lattice QCD configu- 
rations generated by the PACS-CS Collaboration [9] on a 32 3 x 64 lattice using the non- 
perturbatively 0(a)-improved Wilson quark action with Cg^r = 1.715 [Toj and the Iwasaki 
gauge action [ll[ at /3 = 1.90. The aggregate of 2000 MD time units were generated at the 
hopping parameter given by (/c° d , /c°) = (0.13778500,0.13660000), and 80 configurations at 
every 25 MD time units were used for measurements. We then reweight those configura- 
tions to the physical point given by (n u d, K s) = (0.13779625,0.13663375). The reweighting 
shifts the masses of tt and K mesons from = 152(6) MeV and itlk = 509(2) MeV to 
= 135(6) MeV and m# = 498(2) MeV, with the cutoff at the physical point estimated 
to be a- 1 = 2.194(10) GeV 

Observables at the physical point are evaluated through the formula 



(0[U](K ud ,K s )) 



(O[U](n ud ,K s )R ud [U}R s [U]) (K0 



( K udAs) 



( K ud' K s) 



(R uA [U}R s [U]} (K o 



( K ud' K s) 



where the reweighting factors are defined as 



(HI) 



Rud[ 
Rs[ 



det 



det 



D 4[U] J 



(II.2) 
(II.3) 



and LL q [U] is the Wilson-clover quark operator with the hopping parameter n q . We refer to 
Ref. [9] for details of our evaluation of the determinant ratio. Our parameters and statistics 
at the physical point are collected in Table |U 

The relativistic heavy quark formalism 0] is designed to reduce cutoff errors of 0({mQa) n ) 
with arbitrary order n to 0(f (mQa)(ahQCD) 2 ) > once all of the parameters in the relativistic 
heavy quark action are determined nonperturbatively, where f{m,Qa) is an analytic function 
around the massless point rriQa = 0. The action is given by 



(II.4) 
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Dx.y $xy 



- K Q 
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(II.5) 



where kq is the hopping parameter for the heavy quark. The parameters r t ,r s ,CB,CE and 
v are adjusted as follows. We are allowed to choose r t = 1, and we employ a one- loop 
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perturbative value for r s [12J. For the clover coefficients cb and ce, we include the non- 
perturbative contribution in the massless limit Cg^ for three flavor dynamical QCD [loj, and 
calculate the heavy quark mass dependent contribution to one-loop order in perturbation 
theory [l2[ according to 

c b ,e = (c B ,E( m Q a ) ~ Cb,e(0)) PT + Cg W - ( n - 6 ) 

The parameter v is determined non-perturbatively to reproduce the relativistic dispersion 
relation for the spin-averaged IS* states of the charmonium. Writing 

E{ff = E^f + cl s \p\\ (II.7) 

for |p| = 0, (2tt/L), a/2(27t/L), and demanding the effective speed of light c e ff to be unity, 
we find v = 1.1450511 with which we have c e g = 1.002(4). It is noted that the remaining 
cutoff errors are a 2 s f {mQa)(akQ C D) , instead of f(mQa)(aAQ C D) 2 , due to the use of one-loop 
perturbative values in part for the parameters of our heavy quark action. 

We tune the heavy quark hopping parameter to reproduce an experimental value of the 
mass for the spin-averaged IS" states of the charmonium, given by 

M(lS) exp = (M Vc + 3M J/4 ,)/A = 3.0678(3) GeV Q. (II.8) 

This leads to K charm = 0.10959947 for which our lattice QCD measurement yields the value 
M(lS) = 3.067(1) (14) GeV, where the first error is statistical, and the second is a system- 
atic from the scale determination. Our parameters for the relativistic heavy quark action 
are summarized in Table HT1 

We use the following standard operators to obtain meson masses, 

M( g {x) = q f {x)Yq g {x), (II.9) 

where /, g are quark flavors and Y = I, j 5 , 7 M , 27^75, i[7 M , 7^]/2. The meson correlators are 
calculated with a point and exponentially smeared sources and a local sink. The smearing 
function is given by \I>(r) = Aexp(-Br) at r ^ and \1>(0) = 1. We set A = 1.2, B = 0.07 
for the ud quark, A = 1.2, B = 0.18 for the strange quark, and A = 1.2, B = 0.55 for 
the charm quark. The number of source points is quadrupled and polarization states are 
averaged to reduce statistical fluctuations. Statistical errors are analyzed by the jackknife 
method with a bin size of 100 MD time units (4 configurations), as in the light quark 
sector j9[. 

We extract meson masses by fitting correlators with a hyperbolic cosine function. For 
charmonium, Fig. [T] shows effective masses, from which we choose the fitting range to be 
[tmin, tmax] = [10,32]. Similarly, Fig. |2] and Fig. |3] represent effective masses for charmed 
mesons and charmed-strange mesons. We employ the fitting range [t m i n ,t max ] = [14, 20] for 
pseudoscalar mesons, and \t m i n) t max ] = [10,20] for the other channels. 

We calculate the decay constant fps of the heavy-light pseudoscalar meson using the 
improved axial vector current A^ v . 

tf PsPfl = (0\A2»\PS(p)), (11.10) 
A l ™ p = ^2r7 q ^2i^Z A4 {g(x)7475<2(x) 

+c+ 4 «9+ (q(x) l5 Q(x)) + c~ M dl (q(x) l5 Q(x))} , (11.11) 
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where \PS) is the pseudoscalar meson state and d ± is the lattice forward and backward 
derivative. For the renormalization factor Z Aa and the improvement coefficients of the axial 
current c Ai and cL, we employ one-loop perturbation theory to evaluate the mass-dependent 



contributions [13J, adding the nonperturbative contributions in the chiral limit by 



c 



a a ~ (4 4 (m Q a)-4 4 (0)) PT + c™\ (11.12) 



Z M = (Z A4 (m Q a) - Z A4 (0)) PT + Zf, (11.13) 



with cY = -0.03876106 [14| and Zf = 0.781(20) [15|. 



The bare quark mass is determined through the axial vector Ward-Takahashi identity, 

/r)\A imp \P^ 

rr, AWI -L- m AWI - <rn A— 4 ' ' (TT 1 A\ 

m f +7Tlg _ mps ———— } C 1114 ; 

where P is the pseudoscalar meson operator. The renormalized quark mass in the MS 
scheme is given by 

mf (jjl) = Z m (n)mf WI . (11.15) 

Similar to the case of Za 4 , the renormalization factor for the quark mass at the renormal- 
ization scale fi, Z m (fi), is nonperturbatively determined at the massless point, 

Z m (jj) = {Z m {m Q a) - Z m (0)) PT (/i) + Z^(fi), (11.16) 



with Z^ p (fi = l/a) = 1.308(35) [l5|. The charm quark mass is then evolved to fi — ^^arm 
using Nf — 3 four- loop beta function [16j. We employ Nf = 3 based on the fact that our 
simulation includes Nf = 2 + 1 dynamical quarks. 



III. CHARMONIUM SPECTRUM AND CHARM QUARK MASS 

Our results for the charmonium spectrum on the physical point are summarized in Fig. H] 
and Table ITTTl Within the error of 0.5-1%, the predicted spectrum is in reasonable agreement 
with experiment. 

Let us consider the 15 states more closely. Since these states are employed to tune the 
charm quark mass, the central issue here is the magnitude of the hyperfine splitting. Our 
result rrij/^p — m Vc = 0.108(1) (0) GeV, where the first error is statistical and the second 
error is systematic from the scale determination, is 7% smaller than the experimental value 
of 0.117 GeV. In Fig. O we compare the present result on Nf = 2 + 1 flavor dynamical 
configurations with previous attempts on Nf = 2 dynamical and quen ched configurations 
using the same heavy quark formalism and the Iwasaki gluon action [17| . We observe a clear 
trend that incorporation of dynamical light quark effects improves the agreement. 

We should note that the continuum extrapolation is to be performed. A naive order 
counting implies that effects of 0(a 2 s f(mQa)(ahQCD)) from the relativistic heavy quark 
action is at a percent level. Another aspect is that dynamical charm quark effects and 
disconnected loop contributions, albeit reported to give a shift of only a few MeV [l8| . are 
not included in the present work. Additional calculations are needed to draw a definite 
conclusion for the hyperfine splitting of the charmonium spectrum. 
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Using Eq. f ]II.15j) . the charm quark mass is obtained as 



™ c han> = <LJ = 1.260(1)(6)(35) GeV, (III.l) 

where the first error is statistical, the second is systematic from the scale determination, 
and the third from uncertainty in the renormalization factor. The systematic error due to 
the heavy quark of 0(a 2 s f{mQa)(akQCD)) is also to be estimated. Figure |6] compares our 
result with a recent Nf = 2 + 1 lattice QCD estimation by the HPQCD Collaboration [l9[ in 
the continuum limit, which uses the HISQ form of the staggered quark action for the heavy 
quark on the MILC dynamical configurations. 



IV. CHARMED MESON AND CHARMED-STRANGE MESON SPECTRUM 

We calculate the charmed meson and charmed-strange meson masses which are stable on 
our lattice with the spatial size of L = 2.88(1) fm and a lattice cutoff of a -1 = 2.194(10) GeV. 
The D* and D* meson decay channels are not open in our lattice setup. D* and D sl meson 
masses are below the DK threshold |1] but above the D s 7r threshold. Their decays, however, 
are prohibited by the isospin symmetry. On the other hand, Dq and D\ meson masses are 
not computed since their decay channels are open, and therefore a calculation involving Dtt 
contributions is needed. 

Our results are summarized in Fig. [7] and in Table HVl and IVl All our values for the heavy- 
light meson quantities are predictions, because the physical charm quark mass has already 
been fixed with the charmonium spectrum. The experimental spectrum are reproduced in 
2a level. The potential model predicts the D* Q meson mass is above the DK threshold [20], 
which deviates from the experiment significantly. But, our result does not indicate such a 
large difference from the experimental value. A similar result is obtained in other lattice QCD 
calculations [2l[. It should be noticed that our calculation does not cover DK scattering 
states yet. DK contamination for D* and D s \ meson masses can be considerably large. 
Further analysis is required to validate our results for D* and D s \ meson spectrum. 

We compare our results for the hyperfine splittings mjj. —mo and viid s -mj) with 
experiments in Fig. [BJ where we also plot our previous results for Nf = 2 and quenched 
QCD [17|]. The deviation from the experimental value is 1.2<r for charmed mesons, and 2.3a 
for charmed-strange mesons. 

V. CHARMED MESON AND CHARMED-STRANGE MESON DECAY CON- 
STANTS AND CKM MATRIX ELEMENTS 

Table IVT1 presents our estimate of the pseudoscalar decay constants for D and D s mesons. 
Figure M shows the experimental values |l| and our decay constants, as well as three recent 
lattice QCD results: HPQCD and UKQCD Collaboration [HI using HISQ heavy quark on 
the MILC staggered dynamical configurations, Fermilab lattice and MILC group [221 using 
the Fermilab heavy quark on the MILC configurations, and ETM Collaboration [23] who 
uses the twisted mass formalism. Our value for fu s is in accordance with experiment, while 
that for fjj is somewhat larger. Comparing four sets of lattice determinations, we observe, 
both for f D and fn s , an agreement between our values and those of the Fermilab group, 
while there seems to be a discrepancy between our values and those by the HPQCD and 
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UKQCD Collaboration and ETM Collaboration, though continuum extrapolation is needed 
on our part. 

We plot the ratio of fo 3 to fn in Fig. [TUJ Uncertainties coming from the renormalization 
factors cancel out, and that of the lattice cutoff to some extent. Our result is slightly smaller, 
but still Nf = 2 + 1 lattice results are mutually consistent within the errors of a few percent. 

A. Estimating the CKM matrix elements 

The standard model relates \V c d\ to the leptonic decay width of the D meson T(D — > lu) 

by 

F(D -+ lu) = Qftn$m D (l - 4-V \Vcdl 2 , (V.l) 
87T \ mr D ) 

where Gp is the Fermi coupling constant, and mi is the lepton mass in the final state. A 
lattice determination of the D meson decay constant fp with the experimental value of 
T(D — > lu) gives \Vd\- \V CS \ can be obtained in the same way. 

We estimate \V c d\ from our D meson mass and decay constant with the CLEO value of 
F(D — > lu) [24]. Up to our heavy quark discretization error of 0(a 2 s f '(mQa)(aAQCD)) , we 
obtain 

| (lattice) = 0.205(6)(1)(5)(9), (V.2) 

where the first error is statistical, the second is systematic due to the scale determination, the 
third is uncertainty of the renormalization factor, and the forth represents the experimental 
error of the leptonic decay width. For comparison, the PDG value given by \V c d\ = 0.230(11) 
[l| is about 10% larger (see Fig. [TTj) . 

Similarly, using the CLEO value of T(D S — y lu) [25], we find 

|U CS | (lattice) = 1.00(1)(1)(3)(3), (V.3) 

as compared to \V CS \ = 1.02(4) from PDG Q. 

For completeness we also record the ratio |V^ S |/|V^| for which the systematic errors are 
partially dropped out. 

^4 (lattice) = 4.87(14) (0) (0) (27) . (V.4) 
red I 

The PDG value is \V cs \/\V cd \ = 4.45(26). 

VI. CONCLUSION 

We have reported our study of the charm quark system in Nf = 2 + 1 dynamical lattice 
QCD. Although carried out at a finite lattice spacing of a -1 = 2.194(10) GeV, our results for 
the spectra of mesons involving charm quarks are consistent with experiment at a percent 
level, and so are those for the decay constants within a few percent accuracy. These results 
indicate that the heavy quark mass correction rrtQa in the charm quark system is under 
control by the relativistic heavy quark formalism of Ref. jjjj]. Of course, the continuum 
extrapolation and further reductions of statistical noises are required to obtain the result 
competitive with other approaches in the literature. 
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From methodological point of view, we have shown that the realistic heavy quark sim- 
ulations with the light dynamical quark masses precisely tuned to the physical values are 
feasible. With the technique of reweighting, configuration generations are needed to be car- 
ried out approximately around the physical point, and a residual fine tuning to reach the 
physical point only requires a much less time consuming evaluation of the quark determinant 
ratios. Combined with the PACS-CS configuration generation at a smaller lattice spacing 
of gT 1 ~ 3 GeV underway, we hope to return to the issue of continuum extrapolation for 
the charm quark system in future. 
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f3 k uc j k s # conf MD time 

1.90 0.13779625 0.13663375 80 2000 

TABLE I: Simulation parameters. MD time is the number of trajectories multiplied by the trajec- 
tory length. 



K charm v r s C B C E 

0.10959947 1.1450511 1.1881607 1.9849139 1.7819512 
TABLE II: Parameters for the relativistic heavy quark action. 



J r operator lattice experiment 
m^JGeV] Q-+ 75 2.986(1)(13) 2.980(1) 
mj/^GeV] \— -a 3.094(1)(14) 3.097(0) 
~^ e0 [GeV] 0++ / 3.444(33)(15) 3.415(0) 
m Xel [GeV] 1++ 7^75 3.506(30)(15) 3.511(0) 
m fee [GeV] 1+- jaj 3.510(42) (15) 3.525(0) 

TABLE III: Charmonium spectrum in GeV units. The first error is statistical, and the second is 
systematic from the scale determination. Experimental data are also listed [l|. 



J p r operator 


lattice 


experiment 


m D [GeY] 0" 75 


1.871(10)(8) 


1.865(0) 


m D *[GeV] 1" 7i 


1.994(H)(9) 


2.007(0) 



TABLE IV: Charmed meson mass spectrum in GeV units. The first error is statistical, and the 
second is systematic from the scale determination. Experimental data are also listed [l|. 





J p 


r operator 


lattice 


experiment 


m D ,[GeV\ 


0- 


75 


1.958(2)(9) 


1.968(0) 


m D| [GeV] 


1" 


li 


2.095(3)(10) 


2.112(1) 


m D *JGe\] 


0+ 


I 


2.335(35)(10) 


2.318(1) 


m Dsl [GeV] 


1+ 


7*75 


2.451(28)(11) 


2.460(1) 



TABLE V: Charmed-strange meson mass spectrum in GeV units. The first error is statistical, 
and the second is systematic from the scale determination. Experimental data are also listed [l|. 
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lattice experiment 



/p[MeV] 226(6)(1)(5) 206.7(8.9) 
/g a [MeV] 257(2)(1)(5) 257.5(6.1)" 
~W/r> 1.14(3)(0)(0) 1.25(6)~ 

TABLE VI: Our results for decay constants of D meson and D s meson. The first error is statistical, 
the second is systematic from the scale determination, and the third is from the renormalization 
factor. Experimental data are also listed [l|. 
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1.370 



32 x64, Kbd = 0.13779625, k s = 0.13663375 
Kcharm = 0.10959947 




32 x64, Kud = 0.13779625, k$ = 0.13663375 
ICcharm = 0.10959947 





FIG. 1: Effective masses for charmonium. 
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FIG. 3: Effective masses for charmed-strange mesons. 



13 




3) 0.97 

to 



FIG. 4: Our results for the charmonium mass spectrum normalized by the experimental values. 
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FIG. 5: Hyperfine splitting of the charmonium with different number of flavors. 



This \Nork(N/=3 running, a' 1 =2.2 GeV) ^ 
HPQCD(A' / =4 running, a=0) ^ 



1 .24 1^26 



FIG. 6: Comparison of the charm quark mass. The charm quark mass is obtained at Li = a , and 



evolved to fx = fn^ arm using four-loop beta function [la ]. We employ Nf = 3 running based on 



the fact that our simulation includes Nf = 2 + 1 dynamical quarks, while HPQCD collaboration 
uses Nf = 4 reflecting fictitious dynamical charm quark effects 1QQ . 



14 



2 1.03 




2 1.03 



a 1.02 

& 

<n 1.01 
to 

CO 

^ 1.00 



= 0.99 



O 



V, 0.97 
CO 



« 








so D 


n 



FIG. 7: Our results for charmed meson masses(left panel) and charmed-strange meson masses(right 
panel) normalized by the experimental values. 
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FIG. 8: Our results for the hyperfine splittings of charmed meson(left panel) and charmed-strange 
meson(right panel). 
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FIG. 9: Comparison of pseudoscalar decay constants for the charmed meson(left panel) and 
charmed-strange meson (right panel). 
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FIG. 10: Ratios of pseudoscalar decay constants for the charmed meson and charmed-strange 
meson. 
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FIG. 11: Comparison of the CKM matrix elements, |V^|(left panel) and | | (right panel). 
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FIG. 12: Ratio of the CKM matrix elements, \V CS \ and \ V c d 
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